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Abstract 



Recently a new bosonization method has been used to derive, at 
zero fermion density, an effective action for relativistic field theories 
whose partition function is dominated by fermionic composites, chi- 
ral mesons in the case of QCD. This approach shares two important 
features with variational methods: the restriction to the subspace of 
the composites, and the determination of their structure functions by 
a variational calculation. But unlike standard variational methods it 
treats excited states on the same footing as the ground state. 

I extend this method including states of nonvanishing fermion (baryon) 
number and derive an effective action for QCD at finite temperature 
and baryon density. I test the result on a four-fcrmion interaction 
model. 
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1 Introduction 



Increasing temperature and baryon density hadronic matter is expected to 
undergo one or more crossovers and/or phase transitions. Increasing temper- 
ature at zero baryon number one might / should meet a phase in which quarks 
cohexist with hadronic (possibly colored) states pp. Increasing baryon den- 
sity at fixed temperature one should meet a similar phase and possibly a 
color superconducting phase [2] due to a weak attractive channel between 
quarks of different colors. These new states of matter should be at least 
partially accessible to experimental investigation in heavy ion collision ex- 
periments. 

Understanding the behavior of hadronic matter at high temperature and 
baryon density is relevant for the study of early Universe and neutron stars. 
But its theoretical properties can be studied only nonperturbatively and the 
lattice approach, the most powerful tool for first principles, nonperturbative 
studies, is affected in the case of finite density QCD by the well known sign 
problem. 

Some progress was achieved recently O H] by simulations at imaginary 
chemical potential. Other interesting results were obtained [5] within a 
modified version of the Glasgow re- weighting technique and by an approach 
which makes use of a Taylor expansion in the chemical potential in the small 
fi/T region [6]. 

At last a new approach to simulate QCD at finite temperature and 
baryon density was developed [7], which resembles in some aspects that 
of the imaginary chemical potential, but seems to have a wider range of 
applicability [8]. 

The magnitude of quark masses has large effects in numerical simula- 
tions. Several arguments lead to the expectation that the evaluation of 
the fermion determinant is more stable the larger fermion masses are [9]. 
These arguments are relevant to the present work, as explained at the end 
of Section VI. 

To tackle the problem of QCD at finite temperature and baryon density I 
extend a new method constructed to treat fermionic systems whose partition 
function is dominated by fermionic composites. This is certainly the case 
of QCD at low temperature and baryon density, in which the relevant de- 
grees of freedom are mesons and nucleons, and also at high temperature and 
baryon density according to the expectations reported above. This method 
was first developed in the framework of many-body nonrelativistic theories 
[lOj and then applied to relativistic field theories [11] at finite temperature 
and zero fermion density. In the case of QCD, neglecting nucleons it amounts 
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to a bosonization. The heuristic motivation is that reformulation of a the- 
ory in terms of fields related to physical degrees of freedom should make 
it simpler. The starting point of this approach is the partition function in 
operator form, namely the trace of the transfer matrix in the Fock space of 
the fermions. The physical assumption of composite dominance is then im- 
plemented by restricting the trace to fermion composites. This requires an 
approximation of a projection operator on the subspace of the composites, 
the approximation being the better, the higher the number of fermion states 
(called index of nilpotency) in the composites. The approximate projection 
operator is constructed in terms of coherent states of composites, and evalu- 
ation of the trace, which is done exactly, generates a bosonic action in terms 
of the holomorphic variables appearing in the coherent states. 

The structure functions of the composites are determined by a variational 
procedure. So this approach shares two important features with variational 
methods: The restriction to a subspace of the Fock space of fermions, the 
space of chiral mesons in the case of QCD, and the determination of their 
structure functions by a variational calculation. But unlike standard vari- 
ational methods in the present one excited states are treated on the same 
footing as the ground state. 

The utility of variational methods and bosonization has been widely 
appreciated in the theory of many-body systems. But their potentiality 
has also been considered in the framework of relativistic field theories, in 
particular gauge theories, for example by R. Feynman [12] who, however, 
was skeptical about their practical applicability, and recently in connection 
with QCD at high baryon density |13j . 

The approach just outlined is compatible with any regularization. But 
in gauge theories the effective action of the composites will involve vacuum 
expectation values of invariant functions of gauge fields which cannot be 
evaluated within the present framework. Therefore a lattice formulation 
was adopted in order to be able to to extract such expectation values from 
numerical simulations. One is then confronted with the well known difficulty 
with chiral invariance, which can only in part be overcome by using Kogut- 
Susskind fermions. However the method can, at least in principle, be used 
with any other lattice regularization |14j for which a transfer matrix has 
been explicitly constructed. 

The formalism of the transfer matrix does not treat time and space in 
a symmetric way, and therefore Euclidean invariance of the bosonic action 
must be checked a posteriori. All other symmetries are instead respected. 

Before outlining the extension of this approach I will resume what has 
been already done. The validity of the method was tested on a model with a 
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4-fermion interaction in 3+1 dimensions: Euclidean invariance was recovered 
in the continuum limit and all the known results in the boson sector were 
exactly reproduced, namely condensation of a composite boson with the 
right mass, which breaks the discrete chiral invariance of the model. In 
addition the structure function of the composite was determined, and its 
radial factor, in a polar representation, turned out to be identical with that 
of the Cooper pairs of the BCS model of superconductivity. 

To study QCD at finite temperature and baryon density as a first step 
I must introduce quark states in the presence of mesons, namely I must 
construct a Fock space containing composites and their constituents avoiding 
double counting. A similar but more difficult problem has been considered 
since a long time: given a Lagrangian which generates bound states, how to 
replace it by a physically equivalent Lagrangian in which bound states and 
constituents are treated on equal footing [15] . I solve my problem defining 
quasiquark states in such a way that quasiquark-quasiantiquark states are 
orthogonal to meson states. This constraint corresponds to the condition 
on the wave function renormalization of composite particles in the Lehman 
spectral representation of composite operators [15] . 

The next step, the explicit introduction of baryons and antibaryons con- 
structed in terms of quasiquarks and quasiantiquarks is desirable but not 
necessary in a variational calculation, because a space of mesons, quasi- 
quarks and quasiantiquarks obviously contains baryons and antibaryons. 
Therefore in the present paper I will not explicitly include in the parti- 
tion function baryonic states. For a further simplification, which will be 
removed in work in progress, I will exclude antiquarks, so that my varia- 
tional space contains mesons and baryons. This amounts to neglect virtual 
baryons-antibaryons, and it is justified for not too high temperature and 
baryon density. In the resulting effective action the expectation value of the 
chiral sigma field provides a mass to the quasiquarks. On the ground of the 
arguments concerning the effects of quark masses on numerical simulations 
quoted above I hope that numerical simulations with such effective action 
will be more stable. I will also investigate in a separate work the possibility 
of analytical expansions. 

I again test the method on a four-fermion interaction model, reproducing 
the known results in the fermion sector, namely existence of a free fermion 
whose mass is half that of the composite boson, and chiral symmetry restora- 
tion with increasing fermion density. The mechanism of this restoration is 
that quasifermions occupy the lowest energy states, from zero energy up 
to a maximum energy increasing with density, progressively depleting the 
condensate. 
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The paper is organized in the following way. In Section 2 I report the 
general formalism at zero baryon density, in section 3 I define quasiquark 
states and the approximate projection operator in the subspace of mesons 
and quasiquarks, in section 4 I derive a first form of their effective action. 
In Section 5 I apply this action to the study of the four-fermion interaction 
model, deriving the results described above. In section 6 I derive a second 
form of the effective action, which has a more transparent interpretation, and 
allows a crosscheck of the accuracy of the approximation for the projection 
operator by comparison of the results for the four-fermion interaction model, 
which coincide with those obtained by the first effective action. In Section 
7 I summarize my results with an outlook to possible applications. 

2 General formalism at zero baryon density 

To make the paper reasonably selfcontained I report the general formal- 
ism developed for relativistic field theories of fermions in the presence of 
composites dominance at zero fermionic density I make only a small 

modification: I do not fix the gauge, because this would prevent numerical 
simulations and complicate the treatment due to the need to enforce the 
Gauss constraint in Fock space. 

The starting point of this formalism is the standard expression of the 
partition function of QCD in terms of the transfer matrix 

Z = J [dU] exp [S G (U)} Tr F j f[ (f/ V t exp(/xn B )f t+1 ) j (1) 

where Lq is the number of links in the temporal direction, Sq is the gluon 
action and 

Tf = ex p[ — u^Mt it — v^Mjv] exp[-uiVt u\ 

V t = exp[uUnU 0}t u + tfU^ t v}. (2) 

The U^j are matrices whose matrix elements are the link variables at Eu- 
clidean time "t" 

(£W)xi,x 2 = <5 XljX2 C/^ t (xi) . (3) 

Because the formalism treats asymmetrically time and space, I use boldface 
letters, as x, to denote spatial coordinates, and italic letters to denote space- 
time coordinates: x = (i, x). u\ and v\ are, respectively, creation operators 
of quarks and antiquarks in state i, obeying canonical anti-commutation 
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relations. TV is the trace over the Fock space of quarks, /i is the chemical 
potential and ub the baryon number operator. The matrices Mt {Mj being 
the transposed of Mt) and Nt are functions of the spatial link variables at 
time t. They depend on the regularization adopted for the fermions, but 
what follows is not affected by their explicit expressions, which are reported 
in Appendix [B] for Wilson and Kogut-Susskind fermions in the flavor basis. 
I include in the gluon action the term 

5S G = Y,- 4tv - M t ( 4 ) 

t 

which comes from transformations on the fermion fields going from the func- 
tional form to this operator form of the transfer matrix. I introduced the 
notation, which I will use for any matrix A 

tr±A = tr (P (±) A) . (5) 

The operators Pq , which project on the quark antiquark components of 
the quark field are defined in Appendix [Bl tr± is the trace over quarks or 
antiquarks intrinsic quantum numbers and spatial coordinates (but not over 
time). 

The expression (TTJ) for the partition function was given by Luscher [16] 
in the gauge Uq ~ 1, in which Vt = 1 (but one has to impose the Gauss 
constraint in the Fock space of fermions). 

Under the assumption that at low energy the partition function is domi- 
nated by chiral mesons, the trace in the Fock space can be restricted to them. 
The restricted partition function can be written introducing an operator P m 
which projects in the subspace of mesons 

r (No-i ^ 

-2-mesons = 

/ dUexp[-S G (U)]Tr F I f[ (^m 2* Vt exp(/x n B )T m ) \ . 

(6) 

To construct this projector, meson creation operators are introduced 

*U = ^4/^ = E , (7) 

ij 

where x represents their spatial coordinate, K their quantum numbers, like 
radial excitations, spin, flavor, etc, and their structure functions {wave 
functions). 
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Since fermion creation operators are nilpotent, composite creation oper- 
ators & K can be classified according to their index of nilpotency, which is 
the highest integer exponent VL such that 



(8) 



ft counts the number of fermion states in the composite. By analogy with 
systems of elementary bosons coherent states of mesons can be constructed 



exp 



,X,Jf 



|0>, 



(9) 



where the </> x k's are holomorphic variables. But since composites operators 
do not obey canonical commutation relations, the properties of their coher- 
ent states can differ from those of canonical bosonic coherent states. For 
instance the basic property of coherent states cannot be exactly satisfied 



$xK|</>) / </>xK| 



(10) 



However, if the index of nilpotency of the composites is large enough, the 
composites system resembles a canonical bosonic system, and the properties 
of canonical boson coherent states will approximately hold for the composite 
coherent states, as shown in detail in Refs. [10} lllj. 

Hence, under the assumption that the composite operators which domi- 
nate the partition function have a large index of nilpotency, an approximate 
projection operator in the Fock space of the fermions can be defined 



m 





~d(f>d<t>*~ 


/ 


2ni 



where 



2ni 



n 



2ni 



fill 



(12) 



It is important to observe that the space selected by this operator includes 
2 physically equivalent states obtained for (j) = 0, do. They correspond to a 
completely empty or filled lattice. 

The scalar product of coherent states appearing in the definition of the 
projection operator is 



detj 



(13) 
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where 



(14) 



and for any matrix A 

det±A = det(P ± A) . (15) 

X is the identity in the space of all the matrices. I remind that the entries 
of these matrices do not include time. By a little abuse of notation I will 
often write " 1" instead of X. 

With periodic boundary conditions for the gauge fields, the partition 
function is 

Z = Tr F jrj Vo expipns)^ f\ V\ exp(fin B ) ■ ■ ■ ^1 _ 1 V Lo -i exp(^n B )f | 

(16) 

while its restriction to mesons is 



--mesons 



Tr F {p m f$ V t x V m f J V x ■ ■ ■ P m V Lo ^f } 



Lo-l 

n 

t=0 



2iri 



t \f}V t T t+x \ 9t+1 ) 



(9t\9t 



(17) 



where a copy of the Fock space of the mesons has been introduced at each 
time slice. The chemical potential has disappeared because it is not active 
in a space of only mesons. Explicitly 



(18) 



.x.K 



I remark that the structure functions ^x.A' do not depend explicitly on time, 
but as they are functions of gauge fields, time will enter as a label of these 
fields. 

In the evaluation of the trace on the fermionic Fock space the only dif- 
ference with respect to [UJ is the presence of the operator Vt. But I notice 
that the product of operators Vt Tt+x has an expression similar to that of 
f t+ i 

V t f t+1 = exp[-tf ln(e Mt+1 U^ t ) u - v ] ln(e M '+i U^v] exp[vN t u] . (19) 

Then evaluation of the trace over the Fock space proceeds exactly as in [11] 
with the result 



Zmesons= [dU]ex V [S G (U)} / ] [ 



2m 



exp [— iSmesons( 



(20) 



S 



where 



and 



n 



<Pta<Pt 
2iri 



n 

x,K 



dcj)K(x)d(j)* K (x 



2ni 



Smesons — 

t 

In the last equation 



tr_ 



lnR t + ]a72* + M t ' 



(21) 
(22) 



Rt 



(1 + 7*N) t+ i e Mt+1 C/ f it e M * (1 + iV^); 



oM t+ i rrf 



in which I set 



(23) 



(24) 



Notice that the matrix TZt involves gauge fields at time t and t + 1. The 
notation is somewhat different from that of 

It is remarkable that S'mesons has been evaluated exactly [TT], so that 
the only approximations in the partition function are the physical assumption 
of boson dominance and the form of the projector over the meson subspace. 
Since the projector depends on the structure functions $ X ,K) the effective 
action is a functional of these functions which are determined by a variational 
calculation on the quantities of interest. In simple cases, like the four- 
fermion interaction model, the variational calculation provides the exact 
form of the structure function. In QCD, unless some analytic progress is 
made along a way similar to that of the four-fermion interaction model, one 
has to adopt a trial expression. 

In Ref. [TT] an alternative, equivalent form of the effective action was 
derived, which has a more transparent interpretation. I will not report it 
here because I will also derive two forms of the effective action at finite 
baryon density, but for the second one I will follow a somewhat different 
procedure. 



3 Quasiquarks and generalized Bogoliubov trans- 
formations 

In order to extend the formalism to QCD at finite baryon density, I must 
introduce in the partition function states with nonvanishing baryon num- 
ber. In the spirit of composites dominance, I should then construct baryonic 
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composites and define a projection operator on the subspace of mesons and 
baryons. But in the present paper I will be satisfied by introducing, in addi- 
tion to mesons, quasiparticles states with quark quantum numbers, which I 
will call quasiquarks (and quasiantiquarks) . Such a space obviously contains 
the space of mesons and baryons. 

To avoid double counting, mesons and quasiparticles must satisfy a mu- 
tual compositeness condition: Mesonic states must be orthogonal to quasiquark- 
quasiantiquark states. This constraint has the physical meaning of the con- 
dition Z = for bound states in the Lehmann spectral representation of 
composite operators |15j . namely the the condition required to introduce a 
bound state on the same footing as the constituents in a Lagrangian. 

I will denote by quasiparticles destruction operators. I will enforce 
the compositeness condition by requiring that quasiquark-quasiantiquark 
states annihilate coherent states of mesons 



a; 



0. 



(25) 



This condition will prove crucial in the application to the 4-fermion model. 

The compositeness condition in the above form can be solved exactly, 
even though composites coherent states do not satisfy exactly the basic 
property (|10p of coherent states of elementary bosons. The operators Q!j,/3j 
are obtained by a generalized Bogoliubov transformation [T7] 



R5 (u - 
v + tff*) (R' 



(26) 



where 



r' = (i + rj*y 



(27) 



Bogoliubov introduced his tranformation to construct a theory of super- 
conductivity in the presence of an electron-phonon interaction, and defines 
quasiparticles in terms of electron particle-holes states, explicitly breaking 
the f7(l) symmetry related to fermion conservation. The above transforma- 
tions instead respect all symmetries, which is necessary with gauge inter- 
actions, thanks to the presence of the bosonic field 4>. The operators 
and their Hermitean conjugates 



0l 



at 



vF] R2 



(R')s (V 



(28) 
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satisfy canonical commutation relations. As anticipated in the Introduction 
I include only mesons and quasiquarks in my variational space, but not 
quasiantiquarks. My variational space does not contain antibaryons, which 
are not expected to be important for not too high temperature and baryon 
density. In any case it will be clear that extension of the formalism to include 
quasiantiquarks should not present any significant difficulty. 
I define "coherent" states of quasiquarks and mesons 

\a, 4>) = exp(-a • a f ) exp(0 • $ f )|0) (29) 

where the cti are Grassmann variables and 

a ■ a = a i ®i ■ (30) 

i 

These states can be recast in the form 

\a,fa = exp(n t R~5 Q; + . $t)| ) (31) 

which is more convenient for calculations. The operator which approxi- 
mately projects on states of mesons and quasiquarks is 

(a, <p\a, \a, 4>){a, cf>\ (32) 

where the measure is 

{a, (j)\a, (j))" 1 = (a|a:) _1 (0|0) _1 = exp{ tr_ InR — a* • a} . (33) 
If Vm-q is an approximate projector it must satisfy the equations 
(0|6 mi a ni P m _ g (at)"2($tyn 2 | ) „ 

<0| ^^(a^^^lO) cc 5 mi , m2 5 ni , n2 . (34) 
These equations are generated by the following ones 

(4>iai\P m - q \4>2a2) — (0iai|^2a2) , (35) 

by taking derivatives with respect to the variables ai, fa and setting them 
equal to zero. The left hand side of (|35p is 

(cui, 4>i\V m - q \a 2 , 02 ) = J [da*da] 

+ ln(l + J*T X ) + ln(l + t\T) 

+a\ (1 + ^i)- 1 ^ + T\T)- x R" 5 a 2 \ . (36) 



m—q 



[da* da] 



2vri 



2vri 



exp {tr_ [ In R 
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This shows by inspection that the first member of (|34j) vanishes unless mi = 
1^2 , n i = ~n>2- Evaluating the integral in the above equation by the saddle 
point method, as done in [11] for V m , we see that V m - q is approximately a 
projector if we assume 

tr ($t$V ~ Q-n+i _ (37) 
I remind that f2 is the index of nilpotency of $. 



4 First form of the effective action at finite baryon 
density 

I follow the derivation of the effective action outlined in Section 2 for zero 
baryon density. I skip many intermediate steps because calculations of this 
kind have been reported in any detail in |llj . and can be easily repeated here 
by the help of the formulae collected in Appendix [Aj I start by evaluating 
the matrix elements of the transfer matrix between coherent states according 
to 

(a t ,(p t \f^ Vtexp(fih B )f t+1 \at + i,4>t+i) = J [d^dj] [d5*d5] 
xe-^- 5 * s (a t , <f>t\T^8)(j8\Vt exp(^n B )T m |a m , <£ t+1 ) . (38) 
The last factor is 

(j5\ V t exp(fj,n B )T t+ i\a t+ i,(t)t+i) = det_(l + J^N) t +\ 
x exp ^*Uo,te~ Mt+1 (1 + J^N)^ ef R^at+i + 4+i e ~ A/i+1 Kt 5 * } ■ 

(39) 

A similar result for the other matrix element and integration over 7*, 7, 8* , 8 
leads to the expression 



(a t ,(f>t\Tt Vtexp(fj,nB)Tt+i\at+i,(/)t+i) = det_ M *K t 1 
x exp (a\ e» TZ t U Q>t e~ Mt+1 R t 7i «t+i) • ( 40 ) 



From the measure appearing in the definition of V m -q, Eq. (|33j) . I get the 
factor 

(at, <t>t\a-t, (pt)^ 1 = det_ Rt exp (-a* • a t ) . (41) 
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Putting these pieces together I get the effective action of mesons interacting 
with quasiquarks 



S. 



mesons-quarks 



Si 



mesons 



xK,f/„, 1 e- M '+'R 1+ 2 1 a, + i 



-a t + e^Rj 



(42) 



where Smesons is given by Eq. (|22p . I remind that that a is a 2-spinor with 
the quark intrinsic quantum numbers. It can be put in a more transparent 
form 

■^mesons-quarks = Smesons - s ^ a* (V t - H t )a t+1 (43) 

t 

by introducing the lattice covariant derivative in the presence of a chemical 
potential and the lattice Hamiltonian 



V 4 a i+ i 
Ht 



- (e M U ,t att+i ~ a t 
s 

1 



- e 



(44) 



The factor s in the above equations takes the value 2 in the Kogut-Susskind 
regularization because the quarks live on blocks, and 1 in the Wilson regu- 
larization. I notice that the time derivative is not symmetric, so that this 
action does not give rise to fermion doubling. Integrating over the Grass- 
mann variables I get the purely bosonic effective action 



S, 



effective 



Tr_ In [W 1 Vie M 



-.(+) 



(45) 



where I adopted the following notations: All matrices (with exception of 
which do not have a time lable are diagonal in time with matrix 
elements 

(^o)xi,ti,x 2 ,t a = £ii 1 t2*xi,x 8 #0 ) t 1 (xi) 
while the matrix elements of space-time translation operators are 



V / xi,X2 



(47) 



" Tr" is the trace on all entries including time, while I remind that " tr" is 
the trace on intrinsic and spatial quantum numbers only. 
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5 Application to a four-fermion interaction model 
with Kogut-Susskind fermions 

To get insight in the above result and also to test it I apply it to the four- 
fermion interaction model adopted as a test at zero fermion density [TT]. It 
is a model in 3+1 dimensions regularized on a lattice with Kogut-Susskind 
fermions in the flavor basis (I do not know any formulation of the transfer 
matrix in the spin-diagonal basis which can be used in the present formal- 
ism). For each of the four Kogut-Susskind tastes there are Nf degenerated 
flavors. Hence, the continuum limit will describe a theory with 4Nf flavors. 
In the flavor basis the action reads 

s = E ' E t m 1 ® 1 + Qkv Wv) + '$(*M*)) 2 (48) 

x y fx 

where m is the mass parameter, g 2 the coupling constant, ifi the fermion 
fields and Q the hopping matrix: 

Q = ^ 7m ® 11 [pH V W + ?W V(-)] . (49) 

The matrices to the left (right) of the symbol (g) act on Dirac (taste) indices. 
I denote by 7 and t the matrices acting on these indices, respectively. The 
operators 

PjP = ~ [1 ® 1 ± 7^75 ® t 5 ^] • (50) 

are orthogonal projectors. The fermion fields are defined on blocks (see 
Appendix [B] for details). The right and left derivatives are given by 

= 4 C3 ±} - • ( 51 ) 

The factor 1/2 is due to the fact that the operators translate by one 
block. The model has a discreet chiral symmetry at m = 0: 

if) -»• -75 ® *5 ^ : 1p ^ tp"/ 5 ®t 5 . (52) 

To have an action bilinear in the fermion fields a scalar field cr(x) is intro- 
duced, whose integration generates the four-fermion coupling: 

s ' = E ' E t( m + fj ) 1 ® 11 + qumv) + S E ' °*( x ) • ( 53 ) 
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The partition function now reads 



[d(r}[di/;dip] exp [-S'] . 



(54) 



Its restriction to fermion composites plus a fermion gas, without antifermions 
is 



C-F 



J [da][d(ft* d(ft][da* da] exp 



4N 



X 



S c - 



(55) 



Sc-F is given by Eq, (|42j) . in which one has to insert the expressions of 
the matrices M, N appropriate to Kogut-Susskind fermions in the flavor 
basis [18J: The matrix M is equal to zero and the matrix N is reported in 
Appendix [Bj Integration over the fermion fields gives the effective action 



S, 



effective 



Tr_ In 



TT l R — e M T (+) 



(56) 



Now I look for constant values of the fields (ft* , (ft and a which make the 
action stationary. I put a bar over constant fields and their functions. Then 



'effective 



Tr_ In 



TZ R — e^T, 



(+) 



(57) 



and I can perform the sum over time getting 
^effective = ~^ L otr_ j/x. 



K 1 R 



K 1 R-e 



(58) 

where 6 is the step function. For fj, = I recover the effective action de- 
rived [TT] at zero fermion density. 

To determine the magnitude of the condensate I must perform a variation 
with respect to the boson fields (j) , (ft, and to determine the form factor of 
the composite a variation with respect to the matrices , T . But ^effective 
does not depend on these variables separately, it is a function of and T . 
The saddle point equations with respect to these matrices for e M < 1Z R, 
are identical to the ones for zero chemical potential. Using the result of [11] 
I then get 

^ = 4^1 + ^ + 1), e^<n^R, (59) 



2H 



where 



(60) 
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with A given by Eq. (|96p . I notice that H differs from the lattice Hamiltonian 
defined above 



U = ie^ (l - R 1 = e^H (yi + H 2 -H^j , (61) 

but they are equal in the formal limit of vanishing lattice spacing. In this 
limit I can rewrite Sq-f m the form 

S c -F =S c -2j2<**t - {H - M )] 9(2H - n)a t . (62) 

t 

For m = fx = 0, I recover the well known result that the fermionic system 
under consideration in the limit of Nf — > 00 contains free fermions of mass a 
in addition to free bosons of mass 2 a. I emphasize that the result recovered 
in this way is only formal. Indeed after adding one fermion I should deter- 
mine the new minimum of the action, namely the variation of the structure 
functions. But it can be justified in a concrete way by evaluating the dif- 
ference of ^effective given by Eq|64] below at fermion numbers differing by 
one unit. 

Now I impose the condition on the fermion number which determines 
the chemical potential. From Eq. (|58p I get 

2d— r 1 —1 

— — — 5 e |j ec ^ ve = tr_ 6 exp fi — 1Z R = tr_ [exp [i — 1 — 2H] = np . 

(63) 

For [i < 2 a, np = 0. For /j > 2<j, quasifermions occupy the states from zero 
energy up to a maximum energy E np depending on the fermion number np. 
The effective action at the minimum takes the form 



^effective = ~ L o tr - { ln (V 1 + # 2 + H) 2 9 (2 H + 1 - exp fx) j 



(64) 



Stationarity with respect to a yields the gap equation which determines the 
masses and therefore the breaking of chiral invariance 

o = "^effective = 2L o ^ tr - { H ^rjp 9 ( 2 ^ + 1 " exp } ' 

(65) 

Increasing the fermion density, namely the chemical potential, quasifermions 
occupy higher and higher energy states depleting the condensate, until only 
the solution W = remains and chiral invariance is restored. 
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6 Second form of the effective action 



The expression (J59J) of the form factors is somewhat surprising, because 
they are increasing functions of momentum. In [llj a more natural form 
was deduced by performing a unitary transformation in the fermionic Fock 
space and deriving the corresponding effective action. This transformation 
changes the empty lattice into the fully occupied one and particles into holes. 
In this new Fock space the structure functions are decreasing functions of 
momentum, and in a polar representation their polar factor is equal to that 
of the Cooper pairs of the BCS model of superconductivity. 

But in addition the second form of the action provided a test of consis- 
tency of the approximation for the projection operator V m . I could follow 
the same path at nonzero baryon density, but instead I will get a similar 
result in a different way. First I rearrange the trace in Fock space in the 
following way 

Z = Tt f jVo exp(nn B )f 1 T-J V\ exp(/xn B ) • • • Vl -1 ex.p{l^n B )f f^ j . 

(66) 

Then I insert the projection operator "Pm-q hr the trace according to 
^mesons-quarks = f [dU]e^[S G (U)]Tr F \ f[ 



V m -qexp((j,h B )VtTt + ifl +1 ) j . (67) 

I emphasize that Z', Z need not coincide with each other because T-m-q is 
not an exact projection operator, but the results obtained by the two forms 
should agree within the approximation for "Pm-q- A comparison between 
these results provides a check of its accuracy. 

In the same way as for the first form of the effective action I evaluate 
the matrix elements 

{a t , 4> t \ exp(fin B )Vt f t +i f} +1 \a t+ i, <j) t+ \) = exp j-tr+7^ 

+ a* t e» R; 1 U , t e" Mt + 1 K' t+1 a t+1 ) , (68) 



where 



1 + ( N t + e- Mt Ul t _ x T t -x U ,t-i e~ M ' 



t 



x (jv t + e -Mt jr t e -M t ^ -1 e -M} (6Q) 
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Including the contribution (|33[) from the measure I get 



^mesons-quarks = mesons - s £ , 

t 



a* t (V t -H t ) a t+ i 



(70) 



where 



•"mesons 



£ 

t 

1 

-e 

s 



tr_ [-lnR t + ]nK t + Ml 



I'- 



-M t + l -I-,' T> 2 



Integrating over a* , a I get the purely bosonic action 



effective 



-Tr_ In 



-R {K)- l e M + e IM U Q T ( 



(+) 



By exploiting the cyclic property of the trace it can be rewritten 



Si 



effective 



-Tr_ln 



(K ) _i R e 



(71) 



(72) 



(73) 



This expression differs from Effective > Eq. ([^5]) , by the replacement of TZ by 
K'. 

I use this second form of the effective action for the four-fermion interac- 
tion model with Kogut-Susskind fermions in the saddle point approximation. 
By means of the results of ref.pT] I find 



1 



(74) 



Now the structure function is a decreasing function of the constituent fermions 
energy. Using the above expression I find that 



(75) 



so that the results concerning mass of the uncorrelated fermions and restora- 
tion of chiral symmetry derived by the first form of the action are recovered. 

I conclude this Section by an observation about the way the arguments of 
Ref. [9] concerning the stability of numerical simulations might apply to the 
present effective action, comparing QCD with the four-fermion model. In 
this model after linearization there is a field, the sigma field, which appears 
in the matrix N in the same position as the fermion mass, so that its expec- 
tation value provides in a natural way a mass to the fermion in the broken 
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phase. In QCD there is no such field, but the chiral sigma field can play the 
same role. Indeed its expectation value appears in Eq. (|69p as an addendum 
to the matrix N and therefore as an addendum to the bare quark mass. To 
the extent that high quark masses can stabilize numerical simulations, use 
of the present effective action should make these calculations easier. 

7 Summary and outlook 

I extended the formalism of composite boson dominance to the case of non- 
vanishing fermion number. This required a definition of fermion and an- 
tifermion states in the presence of bosonic composites satisfying a composite- 
ness condition to avoid double counting. These fermion (antifermion) states 
are called quasifermions (quasiantifermions). Their definition is achieved by 
a generalized Bogoliubov transformation. 

In the application to QCD I restricted myself to a space of mesons 
and quasiquarks, excluding quasiantiquarks, which contains the space of 
mesons and baryons. Neglecting quasiantiquarks amounts to neglect virtual 
baryons-antibaryons, which is justified for not too high temperature and 
baryon density. Obviously if one wants to investigate any high temperature 
and/or baryon density phase transition quasiantiquarks must be included, 
but I cannot foresee any obstruction in this extension. 

I derived two forms of the effective action and I applied both of them to a 
four-fermi interaction model at zero temperature but finite fermion density. I 
recovered in both ways the known results, which provides a crosscheck of the 
approximation of the projection operator introduced to restrict the fermionic 
Fock space in the partition function. The discrete chiral invariance of the 
model (at zero fermion mass) is broken by composite boson condensation 
and the spectrum of the broken phase contains, in addition to a composite 
boson, a free fermion whose mass is half that of the boson. Increasing the 
fermion density, quasifermions occupy the lowest energy states up to an 
energy which increases with increasing density depleting the condensate, 
until chiral symmetry is restored. The compositeness condition is crucial to 
get these results. Since the action of the composite boson is known [11], one 
could study the system also at finite temperature and density. 

Possible applications of the present formalism include numerical studies 
of the evolution of the state of baryon matter with temperature and density 
and the associated phase transitions. In this connection I remind the way the 
arguments of Ref. [9] concerning the stability of numerical simulations might 
apply: the expectation value of the chiral sigma field appears in Eq. ([69|) as 
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an addendum to the matrix N and therefore as an addendum to the bare 
quark mass. To the extent that high quark masses can stabilize numerical 
simulations, use of the present effective action should make these calculations 
easier. 

Also exotic states of baryon matter can be explored. For instance it is 
not difficult, as it will be shown in a separate paper, to introduce in the 
present formalism diquark states. Among abnormal states of hadronic mat- 
ter I would like to mention the layered spin-isospin phase [21] . This is a state 
with one-dimensional crystallization (which distinguishes it from usual pion 
condensation) in which layers of spin-up protons and spin-down neutrons 
alternate with layers of spin-up neutrons and spin-down protons. An inves- 
tigation of a dynamical realization of such a phase in light deformed nuclei 
showed that the spin-isospin nucleon-nucleon interaction is not sufficiently 
strong to produce it [22] , while the critical density for a static phase in neu- 
tron stars has been estimated [231 to be 3 ~ 4 times normal nuclear density. 
A first principles calculation might nevertheless be worth while to make an 
assessement of some simplifications done in the quoted works. 

To perform numerical simulations it is necessary to adopt a trial expres- 
sion of the mesons structure functions, which should be a function of gauge 
fields depending on temperature and baryon density, as suggested by the ex- 
ample of the four-fermion model. To this end any analytical investigation of 
the effective action, for instance according to an expansion in inverse powers 
of the index of nilpotency, might be of great help. In its absence, the form 
of trial structure functions can be suggested by existing results about the 
spatial structure of hadrons, of which a few examples can be found in [23]. 

Acknowledgments I am grateful to S. Caracciolo and G. DiCarlo for many 
fruitful discussions. This work has been partially supported by EEC under 
the contract "Forces Universe" MRTN-CT-2004- 005104 



A Grassmann integrals and coherent states 

If | a) is a fermionic coherent state 




(76) 



then 




(77) 
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and the identity can be written 

J[dada*] (a\a)- l \a)(a\ = 1. (78) 

I remind the fundamental property of coherent states 

u\a) = a\a) (79) 

which implies the relations 

(ap\exp(vNu)\-f5) = exp{5N"/)(al3\j5} = exp{5N~f + a* 7 + p*5\80) 
(7<5|exp(u t ^" t {) t )|0) = (0|exp(&Fu)|7<f)* = exp( 7 *^ t 5*) . (81) 

With the help of these formulae one can compute matrix elements of the 
type 

(a/3\e em e &i ^\0) (82) 



d-f*d-fd5*d5 



(a/3|e CAr "| 7 <5)(75|e" + ^ t " t |0) (83) 



(7<5|7<5) 

[d 1 *d 1 d5*d5]e-^^- & * &+5N ^ +a ^< + ^ &+ ^ :FU * (84) 

= J [d5*d5] e-W+r^t+FS-r™ (85) 
= exp {tr_ ln(l + F*N T ) - (3*(1 + F*N T )- l F*a*} , (86) 
by use of the identity 

J [da* da] exp(-aMa + J* a + a* J) = del A exp(JM _1 J) . (87) 

B The matrices M, N of the transfer matrix 

In this Appendix I report the expressions of the matrices M, N appearing in 
the definition of the transfer matrix for the Kogut-Susskind and the Wilson 
regularization. Their common feature is that they depend only on the spatial 
link variables. 

B.l Kogut-Susskind's regularization 

Kogut-Susskind fermions in the flavor basis are defined on hypercubes whose 
sides are twice the basic lattice spacing. While in the text intrinsic quantum 
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numbers and spatial coordinates were comprehensively represented by one 
index i, here I distinguish the spinorial index a = {1, . . . , 4}, the taste index 
a = {1, . . . , 4} and the flavour index i={l, Nf}, while x = {t, xi, . . . , x^} 
is a 4- vector of even integer coordinates ranging in the intervals [0, Lt — 1] 
for the time component and [0, L s — 1] for each of the spatial components. 
I distinguish summations over basic lattice and hypercubes according to 

X X 

The projection operators over fermions-antifermion states are 

P (±) = ^(1 ® 1 ± 7075 ® t 5 to) ■ (89) 

The relation between the variables u, v and the quark q field is 

Pi ) +) q=\u, P^q = -^. (90) 

In the presence of the scalar field a and of gauge fields, neglecting an irrel- 
evant constant, M = 0, while iV is 1181 



P H v (+) + p(+) v H 

J J ^ 3 3 



N = -2 ^ (m + o-)7o ® 1 + E 7o7i ® 1 



where 

V< +) = \{U 3 T^-1) (91) 

V« = ^(l-^-^j) 02) 

are the lattice covariant derivative as the are the forward and back- 
ward translation operators of one block, that is of two lattice spacing in the 
original lattice, in the ji direction (with unit versor ft) 



I set 



P/i ]a;i,a:2 ~~ ^x 2 ,x 1 ±2jl ■ (93) 

N ] N = AH 2 . (94) 



In the absence of gauge fields 

H 2 = (m + a) 2 - A (95) 
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with 

A 

! 

i=l,3 

The eigenvalues of H 2 are therefore the fermion energies 



1 £ (96) 



E 2 = m 2 +p 2 , (97) 
where momentum component pj is 

p 2 i=\(l-cos2 Pl ). (98) 

and 

P 2 = X> 2 (99) 
i=i 

B.2 Wilson's regularization 

The projection operators over fermions-antifermions are 

P (±) = i(l± 70 ). (100) 

in a basis in which 70 = diag(l, 1, —1, —1). 

The relations between the quark field q and its upper and lower compo- 
nents u, v are 

P { Q +) q = B-^u, P (_) g = B-3«t , (101) 

where 

3 

and K is the hopping parameter. The matrices M, N are 

M = -Ilnf — 
N = 2K cB~^ , (103) 

where 

j'=i 
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